3.151 \(\int (a+b \sin (c+d x))^2 \tan (c+d x) \, dx\)
Optimal. Leaf size=78 \[ -\frac{2 a b \sin (c+d x)}{d}-\frac{(a-b)^2 \log (\sin (c+d x)+1)}{2 d}-\frac{(a+b)^2 \log (1-\sin (c+d x))}{2 d}-\frac{b^2 \sin ^2(c+d x)}{2 d} \]
[Out]
-((a + b)^2*Log[1 - Sin[c + d*x]])/(2*d) - ((a - b)^2*Log[1 + Sin[c + d*x]])/(2*d) - (2*a*b*Sin[c + d*x])/d -
(b^2*Sin[c + d*x]^2)/(2*d)
________________________________________________________________________________________
Rubi [A] time = 0.0709252, antiderivative size = 78, normalized size of antiderivative = 1.,
number of steps used = 6, number of rules used = 4, integrand size = 19, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.21, Rules used =
{2721, 801, 633, 31} \[ -\frac{2 a b \sin (c+d x)}{d}-\frac{(a-b)^2 \log (\sin (c+d x)+1)}{2 d}-\frac{(a+b)^2 \log (1-\sin (c+d x))}{2 d}-\frac{b^2 \sin ^2(c+d x)}{2 d} \]
Antiderivative was successfully verified.
[In]
Int[(a + b*Sin[c + d*x])^2*Tan[c + d*x],x]
[Out]
-((a + b)^2*Log[1 - Sin[c + d*x]])/(2*d) - ((a - b)^2*Log[1 + Sin[c + d*x]])/(2*d) - (2*a*b*Sin[c + d*x])/d -
(b^2*Sin[c + d*x]^2)/(2*d)
Rule 2721
Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*tan[(e_.) + (f_.)*(x_)]^(p_.), x_Symbol] :> Dist[1/f, Subst[I
nt[(x^p*(a + x)^m)/(b^2 - x^2)^((p + 1)/2), x], x, b*Sin[e + f*x]], x] /; FreeQ[{a, b, e, f, m}, x] && NeQ[a^2
- b^2, 0] && IntegerQ[(p + 1)/2]
Rule 801
Int[(((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_)))/((a_) + (c_.)*(x_)^2), x_Symbol] :> Int[ExpandIntegrand[(
(d + e*x)^m*(f + g*x))/(a + c*x^2), x], x] /; FreeQ[{a, c, d, e, f, g}, x] && NeQ[c*d^2 + a*e^2, 0] && Integer
Q[m]
Rule 633
Int[((d_) + (e_.)*(x_))/((a_) + (c_.)*(x_)^2), x_Symbol] :> With[{q = Rt[-(a*c), 2]}, Dist[e/2 + (c*d)/(2*q),
Int[1/(-q + c*x), x], x] + Dist[e/2 - (c*d)/(2*q), Int[1/(q + c*x), x], x]] /; FreeQ[{a, c, d, e}, x] && NiceS
qrtQ[-(a*c)]
Rule 31
Int[((a_) + (b_.)*(x_))^(-1), x_Symbol] :> Simp[Log[RemoveContent[a + b*x, x]]/b, x] /; FreeQ[{a, b}, x]
Rubi steps
\begin{align*} \int (a+b \sin (c+d x))^2 \tan (c+d x) \, dx &=\frac{\operatorname{Subst}\left (\int \frac{x (a+x)^2}{b^2-x^2} \, dx,x,b \sin (c+d x)\right )}{d}\\ &=\frac{\operatorname{Subst}\left (\int \left (-2 a-x+\frac{2 a b^2+\left (a^2+b^2\right ) x}{b^2-x^2}\right ) \, dx,x,b \sin (c+d x)\right )}{d}\\ &=-\frac{2 a b \sin (c+d x)}{d}-\frac{b^2 \sin ^2(c+d x)}{2 d}+\frac{\operatorname{Subst}\left (\int \frac{2 a b^2+\left (a^2+b^2\right ) x}{b^2-x^2} \, dx,x,b \sin (c+d x)\right )}{d}\\ &=-\frac{2 a b \sin (c+d x)}{d}-\frac{b^2 \sin ^2(c+d x)}{2 d}+\frac{(a-b)^2 \operatorname{Subst}\left (\int \frac{1}{-b-x} \, dx,x,b \sin (c+d x)\right )}{2 d}+\frac{(a+b)^2 \operatorname{Subst}\left (\int \frac{1}{b-x} \, dx,x,b \sin (c+d x)\right )}{2 d}\\ &=-\frac{(a+b)^2 \log (1-\sin (c+d x))}{2 d}-\frac{(a-b)^2 \log (1+\sin (c+d x))}{2 d}-\frac{2 a b \sin (c+d x)}{d}-\frac{b^2 \sin ^2(c+d x)}{2 d}\\ \end{align*}
Mathematica [A] time = 0.13574, size = 64, normalized size = 0.82 \[ -\frac{4 a b \sin (c+d x)+(a-b)^2 \log (\sin (c+d x)+1)+(a+b)^2 \log (1-\sin (c+d x))+b^2 \sin ^2(c+d x)}{2 d} \]
Antiderivative was successfully verified.
[In]
Integrate[(a + b*Sin[c + d*x])^2*Tan[c + d*x],x]
[Out]
-((a + b)^2*Log[1 - Sin[c + d*x]] + (a - b)^2*Log[1 + Sin[c + d*x]] + 4*a*b*Sin[c + d*x] + b^2*Sin[c + d*x]^2)
/(2*d)
________________________________________________________________________________________
Maple [A] time = 0.039, size = 82, normalized size = 1.1 \begin{align*} -{\frac{{a}^{2}\ln \left ( \cos \left ( dx+c \right ) \right ) }{d}}-2\,{\frac{ab\sin \left ( dx+c \right ) }{d}}+2\,{\frac{ab\ln \left ( \sec \left ( dx+c \right ) +\tan \left ( dx+c \right ) \right ) }{d}}-{\frac{ \left ( \sin \left ( dx+c \right ) \right ) ^{2}{b}^{2}}{2\,d}}-{\frac{{b}^{2}\ln \left ( \cos \left ( dx+c \right ) \right ) }{d}} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
int((a+b*sin(d*x+c))^2*tan(d*x+c),x)
[Out]
-1/d*a^2*ln(cos(d*x+c))-2*a*b*sin(d*x+c)/d+2/d*a*b*ln(sec(d*x+c)+tan(d*x+c))-1/2*b^2*sin(d*x+c)^2/d-1/d*b^2*ln
(cos(d*x+c))
________________________________________________________________________________________
Maxima [A] time = 1.63094, size = 95, normalized size = 1.22 \begin{align*} -\frac{b^{2} \sin \left (d x + c\right )^{2} + 4 \, a b \sin \left (d x + c\right ) +{\left (a^{2} - 2 \, a b + b^{2}\right )} \log \left (\sin \left (d x + c\right ) + 1\right ) +{\left (a^{2} + 2 \, a b + b^{2}\right )} \log \left (\sin \left (d x + c\right ) - 1\right )}{2 \, d} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate((a+b*sin(d*x+c))^2*tan(d*x+c),x, algorithm="maxima")
[Out]
-1/2*(b^2*sin(d*x + c)^2 + 4*a*b*sin(d*x + c) + (a^2 - 2*a*b + b^2)*log(sin(d*x + c) + 1) + (a^2 + 2*a*b + b^2
)*log(sin(d*x + c) - 1))/d
________________________________________________________________________________________
Fricas [A] time = 1.60878, size = 186, normalized size = 2.38 \begin{align*} \frac{b^{2} \cos \left (d x + c\right )^{2} - 4 \, a b \sin \left (d x + c\right ) -{\left (a^{2} - 2 \, a b + b^{2}\right )} \log \left (\sin \left (d x + c\right ) + 1\right ) -{\left (a^{2} + 2 \, a b + b^{2}\right )} \log \left (-\sin \left (d x + c\right ) + 1\right )}{2 \, d} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate((a+b*sin(d*x+c))^2*tan(d*x+c),x, algorithm="fricas")
[Out]
1/2*(b^2*cos(d*x + c)^2 - 4*a*b*sin(d*x + c) - (a^2 - 2*a*b + b^2)*log(sin(d*x + c) + 1) - (a^2 + 2*a*b + b^2)
*log(-sin(d*x + c) + 1))/d
________________________________________________________________________________________
Sympy [F] time = 0., size = 0, normalized size = 0. \begin{align*} \int \left (a + b \sin{\left (c + d x \right )}\right )^{2} \tan{\left (c + d x \right )}\, dx \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate((a+b*sin(d*x+c))**2*tan(d*x+c),x)
[Out]
Integral((a + b*sin(c + d*x))**2*tan(c + d*x), x)
________________________________________________________________________________________
Giac [B] time = 5.72047, size = 10604, normalized size = 135.95 \begin{align*} \text{result too large to display} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate((a+b*sin(d*x+c))^2*tan(d*x+c),x, algorithm="giac")
[Out]
-1/4*(4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)^4*tan(1/2*c) + 2*tan(1/2*d*
x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2*tan(1/2*d*x)*tan(1/2
*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/
2*c)^2*tan(c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/2*c) -
2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*tan(1/2*
d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(d*x)^2*tan(1/2*d
*x)^2*tan(1/2*c)^2*tan(c)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)
^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(c)^2 + 2*b^2*log
(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*t
an(c) + 1))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(c)^2 - b^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(c
)^2 + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)^4*tan(1/2*c) + 2*tan(1/2*d*
x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2*tan(1/2*d*x)*tan(1/2
*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/
2*c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/2*c) - 2*tan(1/2
*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*tan(1/2*d*x)*tan(
1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(d*x)^2*tan(1/2*d*x)^2*tan
(1/2*c)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(
d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(tan(d*x)
^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan(
1/2*d*x)^2*tan(1/2*c)^2 + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)^4*tan(1
/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2*
tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*tan(d*x)^2*t
an(1/2*d*x)^2*tan(c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/
2*c) - 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*t
an(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(d*x)^2*ta
n(1/2*d*x)^2*tan(c)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan
(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan(1/2*d*x)^2*tan(c)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(
tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x
)^2*tan(1/2*d*x)^2*tan(c)^2 - 16*a*b*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)*tan(c)^2 + 4*a*b*log(2*(tan(1/2*c)^2
+ 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)^4*tan(1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x
)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(
1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*tan(d*x)^2*tan(1/2*c)^2*tan(c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 +
1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/2*c) - 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4
+ 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2
*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(d*x)^2*tan(1/2*c)^2*tan(c)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan
(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2
*tan(1/2*c)^2*tan(c)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*ta
n(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan(1/2*c)^2*tan(c)^2 - 16*a*b*tan(d*x)^2*tan(1/2*d*x
)*tan(1/2*c)^2*tan(c)^2 + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)^4*tan(1
/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2*
tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*tan(1/2*d*x)
^2*tan(1/2*c)^2*tan(c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(
1/2*c) - 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2
*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(1/2*d*x
)^2*tan(1/2*c)^2*tan(c)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2
*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(c)^2 + 2*b^2*log(4*(tan(c)^2
+ 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*t
an(1/2*d*x)^2*tan(1/2*c)^2*tan(c)^2 + b^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 4*b^2*tan(d*x)*tan(1/2*d*x)
^2*tan(1/2*c)^2*tan(c) - b^2*tan(d*x)^2*tan(1/2*d*x)^2*tan(c)^2 - b^2*tan(d*x)^2*tan(1/2*c)^2*tan(c)^2 + b^2*t
an(1/2*d*x)^2*tan(1/2*c)^2*tan(c)^2 + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*
d*x)^4*tan(1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2
*d*x)^3 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*
tan(d*x)^2*tan(1/2*d*x)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan
(1/2*c) - 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 -
2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(d*x)^2
*tan(1/2*d*x)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2
+ tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan(1/2*d*x)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(
c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan(1/2*d*x
)^2 - 16*a*b*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c) + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^
2 + 2*tan(1/2*d*x)^4*tan(1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)
^2 - 2*tan(1/2*d*x)^3 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan
(1/2*c) + 1))*tan(d*x)^2*tan(1/2*c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/
2*d*x)^4*tan(1/2*c) - 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1
/2*d*x)^3 - 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)
)*tan(d*x)^2*tan(1/2*c)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2
*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan(1/2*c)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(tan(d*x
)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan
(1/2*c)^2 - 16*a*b*tan(d*x)^2*tan(1/2*d*x)*tan(1/2*c)^2 + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1
/2*c)^2 + 2*tan(1/2*d*x)^4*tan(1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(
1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) -
2*tan(1/2*c) + 1))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2
- 2*tan(1/2*d*x)^4*tan(1/2*c) - 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2
+ 2*tan(1/2*d*x)^3 - 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1
/2*c) + 1))*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c
) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*b^2*log(4*(tan(
c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) +
1))*tan(1/2*d*x)^2*tan(1/2*c)^2 + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)
^4*tan(1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x
)^3 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*tan(
d*x)^2*tan(c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/2*c) -
2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*tan(1/2*
d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(d*x)^2*tan(c)^2
+ 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2
*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan(c)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*ta
n(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2*tan(c)^2 + 16*a*b*tan(d*x)^2*tan(
1/2*d*x)*tan(c)^2 + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)^4*tan(1/2*c)
+ 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2*tan(1/
2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*tan(1/2*d*x)^2*tan
(c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/2*c) - 2*tan(1/2*
d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*tan(1/2*d*x)*tan(1
/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(1/2*d*x)^2*tan(c)^2 + 2*a^
2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d
*x)*tan(c) + 1))*tan(1/2*d*x)^2*tan(c)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(
c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(1/2*d*x)^2*tan(c)^2 + 16*a*b*tan(d*x)^2*ta
n(1/2*c)*tan(c)^2 - 16*a*b*tan(1/2*d*x)^2*tan(1/2*c)*tan(c)^2 + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4
*tan(1/2*c)^2 + 2*tan(1/2*d*x)^4*tan(1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^
2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*
d*x) - 2*tan(1/2*c) + 1))*tan(1/2*c)^2*tan(c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2
- 2*tan(1/2*d*x)^4*tan(1/2*c) - 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2
+ 2*tan(1/2*d*x)^3 - 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1
/2*c) + 1))*tan(1/2*c)^2*tan(c)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + ta
n(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(1/2*c)^2*tan(c)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(t
an(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(1/2*
c)^2*tan(c)^2 - 16*a*b*tan(1/2*d*x)*tan(1/2*c)^2*tan(c)^2 + b^2*tan(d*x)^2*tan(1/2*d*x)^2 + b^2*tan(d*x)^2*tan
(1/2*c)^2 - b^2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 4*b^2*tan(d*x)*tan(1/2*d*x)^2*tan(c) + 4*b^2*tan(d*x)*tan(1/2*c)
^2*tan(c) - b^2*tan(d*x)^2*tan(c)^2 + b^2*tan(1/2*d*x)^2*tan(c)^2 + b^2*tan(1/2*c)^2*tan(c)^2 + 4*a*b*log(2*(t
an(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)^4*tan(1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 +
tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*
x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*tan(d*x)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2
*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/2*c) - 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/
2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*t
an(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(d*x)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*t
an(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(ta
n(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(d*x)^
2 + 16*a*b*tan(d*x)^2*tan(1/2*d*x) + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d
*x)^4*tan(1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*
d*x)^3 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*t
an(1/2*d*x)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/2*c) - 2*
tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*tan(1/2*d*
x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(1/2*d*x)^2 + 2*a^2
*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*
x)*tan(c) + 1))*tan(1/2*d*x)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d
*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(1/2*d*x)^2 + 16*a*b*tan(d*x)^2*tan(1/2*c) - 16*a*b*t
an(1/2*d*x)^2*tan(1/2*c) + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)^4*tan(
1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2
*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*tan(1/2*c)^
2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/2*c) - 2*tan(1/2*d*x)
^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*tan(1/2*d*x)*tan(1/2*c
)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1))*tan(1/2*c)^2 + 2*a^2*log(4*(tan(c)
^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1)
)*tan(1/2*c)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 +
tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(1/2*c)^2 - 16*a*b*tan(1/2*d*x)*tan(1/2*c)^2 + 4*a*b*log(2*(tan(1/2*c
)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)^4*tan(1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*
d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^3 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + t
an(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1))*tan(c)^2 - 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*ta
n(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/2*c) - 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*t
an(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x
) + 2*tan(1/2*c) + 1))*tan(c)^2 + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(
d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(c)^2 + 2*b^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4*tan(c
)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1))*tan(c)^2 + 16*a*b*tan(1
/2*d*x)*tan(c)^2 + 16*a*b*tan(1/2*c)*tan(c)^2 + b^2*tan(d*x)^2 - b^2*tan(1/2*d*x)^2 - b^2*tan(1/2*c)^2 + 4*b^2
*tan(d*x)*tan(c) + b^2*tan(c)^2 + 4*a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 + 2*tan(1/2*d*x)
^4*tan(1/2*c) + 2*tan(1/2*d*x)^3*tan(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x
)^3 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)) - 4*
a*b*log(2*(tan(1/2*c)^2 + 1)/(tan(1/2*d*x)^4*tan(1/2*c)^2 - 2*tan(1/2*d*x)^4*tan(1/2*c) - 2*tan(1/2*d*x)^3*tan
(1/2*c)^2 + tan(1/2*d*x)^4 + 2*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^3 - 2*tan(1/2*d*x)*tan(1/2*c)^2 +
2*tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)) + 2*a^2*log(4*(tan(c)^2 + 1)/(tan(d*x)^4
*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c) + 1)) + 2*b^2*log(4*(ta
n(c)^2 + 1)/(tan(d*x)^4*tan(c)^2 - 2*tan(d*x)^3*tan(c) + tan(d*x)^2*tan(c)^2 + tan(d*x)^2 - 2*tan(d*x)*tan(c)
+ 1)) + 16*a*b*tan(1/2*d*x) + 16*a*b*tan(1/2*c) - b^2)/(d*tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(c)^2 + d*
tan(d*x)^2*tan(1/2*d*x)^2*tan(1/2*c)^2 + d*tan(d*x)^2*tan(1/2*d*x)^2*tan(c)^2 + d*tan(d*x)^2*tan(1/2*c)^2*tan(
c)^2 + d*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(c)^2 + d*tan(d*x)^2*tan(1/2*d*x)^2 + d*tan(d*x)^2*tan(1/2*c)^2 + d*ta
n(1/2*d*x)^2*tan(1/2*c)^2 + d*tan(d*x)^2*tan(c)^2 + d*tan(1/2*d*x)^2*tan(c)^2 + d*tan(1/2*c)^2*tan(c)^2 + d*ta
n(d*x)^2 + d*tan(1/2*d*x)^2 + d*tan(1/2*c)^2 + d*tan(c)^2 + d)